Energy-Enstrophy Stability of /3-plane Kolmogorov Flow with Drag 
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We develop a new nonlinear stability method, the Energy-Enstrophy (EZ) method, that is specialized to two- 
dimensional hydrodynamics; the method is applied to a /3-plane flow driven by a sinusoidal body force, and 
retarded by drag with damping time-scale p . The standard energy method (Fukuta and Murakami, J. Phys. 
Soc. Japan, 64, 1995, pp 3725) shows that the laminar solution is monotonically and globally stable in a cer- 
tain portion of the (p, /3)-parameter space. The EZ method proves nonlinear stability in a larger portion of 
the (p, /3)-parameter space. And by penalizing high wavenumbers, the EZ method identifies a most strongly 
amplifying disturbance that is more physically realistic than that delivered by the energy method. Linear insta- 
bility calculations are used to determine the region of the (p, /3)-parameter space where the flow is unstable to 
infinitesimal perturbations. There is only a small gap between the linearly unstable region and the nonlinearly 
stable region, and full numerical solutions show only small transient amplification in that gap. 
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I. INTRODUCTION 

Kolmogorov flow is the simplest example of a two- 
dimensional motion forced at a single spatial scale. This pro- 
vides an opportunity to understand the implications of the dual 
conservation laws of energy and enstrophy, the first of which 
is Fj0rtoft's observation that nonlinear interactions transfer 
energy simultaneously up and down scale QUI]. 

The signature of Kolmogorov flow is that motion is main- 
tained against dissipation by a sinusoidal body force. Kol- 
mogorov viewed his problem as an idealized example of a 
forced-dissipative system in which it might be possible to un- 
derstand the sequence of bifurcations resulting from increas- 
ing the Reynolds number [3]. The first steps in this pro- 
gram were an analysis of the viscous linear stability problem 
04], [5|]. The weakly nonlinear theory, pivoted about the critical 
Reynolds number \/2, was subsequently developed^^]. 

Laboratory experiments using either soap films Hi], or 
shallow layers 11011 . can be driven by electromagnetic forcing, 
or by the motion of an enveloping gas, so that the body force 
approximates a sinusoid. Thus Kolmogorov's problem is also 
important as an experimentally accessible flow in which as- 
pects of two-dimensional hydrodynamics can be tested. In the 
laboratory the main dissipative mechanism is drag on the ad- 
jacent walls, rather than lateral viscosity UHO- 

In the geophysical context, the instability of planetary 
waves on a /3-plane M Q [l5ll is similar to Kolmogorov's 
problem in many respects. Ekman friction, which is equiva- 
lent to sidewall drag in the laboratory, also has a strong effect 
on the stability of planetary waves [ 16]. 

Another geophysical motivation for studying the Kol- 
mogorov problem is the equilibration of baroclinic turbulence. 
The most unstable mode of baroclinic instability is an expo- 
nentially growing sinusoidal flow, that is an exact solution of 
the equations of motion. The amplitude of this mode is lim- 
ited by a secondary instability, resembling Kolmogorov in- 
stability. By deflecting energy into the barotropic mode, and 
thence into zonal jets, this secondary instability equilibrates 
baroclinic turbulence by direct cascade of the thermal mode 
to high wavenumbers 01711 . Thus weakly nonlinear /3-plane 
Kolmogorov flow has been studied as a model of zonal jet for- 



mation in the geophysical context f 1811 . 

Our main concern in this investigation is nonlinear stabil- 
ity analysis of Kolmogorov flow. This avenue was opened by 
Fukuta and Murakami 11911 using the energy method 12(1 12 ill . 
The energy method provides a sufficient condition for non- 
linear stability by finding the critical value of the dissipation 
ensuring that the disturbance energy decreases monotonically 
to zero. Our interest in this question is whether the second 
two-dimensional conservation law, namely squared vorticity 
or enstrophy, might be used to improve the nonlinear stability 
results in 11911 . 

In Sec. |II]we formulate the Kolmogorov stability problem. 
In Sec. [nl]we discuss the linear stability of the flow, focus- 
ing on the limit in which the drag is much stronger than vis- 
cosity. In Sec. |IV| we extend the energy-stability condition 
of Fukuta and Murakami [19] to the /3-plane and thus obtain 
a sufficient condition for nonlinear stability. Comparing the 
results of Sec. [Ill] with those of Sec. [IV] we see that there 
is a region of parameter space in which the flow is linearly 
stable, but the energy method fails to prove nonlinear stabil- 
ity. In Sec. [V] we develop a new nonlinear stability method, 
the Energy-Enstrophy (EZ) method, which is specialized to 
two-dimensional hydrodynamics. The EZ method provides a 
sufficient condition for stability which is stronger than the en- 
ergy method, and consequently the gap between the results of 
linear stability and the nonlinearly stable region of parameter 
space is narrowed, but not eliminated. Section fVl\ concludes 
the paper. 



II. FORMULATION OF THE STABILITY PROBLEM 



With non-dimensional variables, the vorticity equation is 



V 2 -0t + JO, V 2 V>) + (3ip x = vV 4 ip - pV 2 ip + cos(x - 



x f) 



(1) 

In Eq. ([TJ the incompressible velocity field is obtained from 
a stream function ip(x,t) according to (it, v) — (—i/j y ,ip x ). 
The domain is a doubly periodic square 2ttL x 2ttL, where 
L is an integer. The relative vorticity is V 2 ^, where V 2 = 
<9 2 + dy is the two-dimensional Laplacian, and J(a. b) = 
the Jacobian. /3 is the gradient of the Cori- 
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olis parameter along y and the two dissipative mechanisms, 
viscosity v and drag \i, are represented by the first and second 
terms on the right hand side of Eq. (Q~|). 

The flow in Eq. (Q} is forced by a sinusoidal body-force, 
which is the signature of the Kolmogorov flow. In dimensional 
variables, the Kolmogorov forcing is specified as 

tJ 2 cos kf(x — Xf) ; (2) 

there is a length scale kj 1 and a time scale r/. To obtain the 
non-dimensional form in Eq. ((TJ, we have scaled using r/ and 
kj. Thus, if * denotes a dimensional quantity, then in Eq. (UJ 
the non-dimensional control parameters are (3 = Tfj3*/kj, 
[i = TffJ,*, v = kfTfi;* and L = fc/L*. Also in Eq. ((T), the 
phase Xf = tan -1 [/3/(/i + v)\ is defined so that there is a 
steady laminar solution: 

V'l(x) = ~ acosx , (3) 

where the amplitude of the laminar flow is 

a{p,n,v)= (4) 
A. Dynamics of the disturbance 




time 



FIG. 1: (Color online) Disturbance energy E v (t) from the numeri- 
cal solution of Eq. (TJ with f3 = and v — 10 -3 . The results are 
normalized by the laminar energy E\_ = \a z . The solid and dashed 
curves illustrate the transient growth in E v when the laminar solu- 
tion ipi_(x) is globally stable, but not monotonically globally stable. 
The dot-dashed curve illustrates monotonic global stability, and the 
dotted curve shows the development of a linearly unstable perturba- 
tion. 



The disturbance (p(x,t) to the laminar stream function is 
defined by 

tp(x,t) = i/)l(x) + (p{x,t) , (5) 

and the equation of motion of cp is obtained by substituting 
Eq. © into Eq. ([JJ 

V 2 tpt + J(ip, V 2 ^) + asinx (V 2 tp + ip) + (3ip x 

= fVV-/jV 2 p. (6) 

The disturbance energy E v and enstrophy are defined 

as 

E V ^\(\^V\ 2 ) and Z v = i((VV) a >, (7) 

where (• • • ) denotes spatial average over the whole domain. 
Multiplying Eq. © by ip and averaging gives the disturbance 
energy equation: 

dE 

= a (<p x <py cos x) — 2fiE<p — 1vZ v . (8) 

Multiplying by V 2 (p and averaging produces the disturbance 
enstrophy equation: 

= a {ipxVy cos x) - 2[iZ v - 2vP v , (9) 

where P v = | (|VV 2 <^| 2 ) is the disturbance palinstrophy. 

It is remarkable that the disturbance energy and disturbance 
enstrophy are generated at equal rates i.e., the same source, 



namely a (f x fy cos x), appears in Eq. (0 and in Eq. (0. Thus 
subtracting Eq. (O from Eq. (0 gives 

~(E V - Z v ) = -lx{E ¥ - Z v ) - vZ v + vP v . (10) 

This cancellation is a general property of flows forced by a 
single Helmholtz eigenmode, and is the basis for recent con- 
straints on the spectral distribution of energy and enstrophy 
in two-dimensional turbulence 11221 l23l 12411 . If v = then 
solution of Eq. (fTOb shows that the difference E v (t) — Z v (t) 
decays exponentially to zero. This observation is the basis of 
the EZ method in section [V] 



B. Linear instability, global stability and monotonic global 
stability 

Following Joseph [20], we say that the laminar solution 
ij}\_{x) is globally stable if for all initial disturbances 

lim£L(i)=0. (11) 

An even stronger form of stability is monotonic global stabil- 
ity, meaning that dE v /dt < for all t > 0. In both cases the 
laminar solution ultimately attracts all initial conditons. 

Global stability does not forbid transient increases in E v (t) 
e.g., see Fig. [JJ It is a limitation of the standard energy 
method, reviewed in section IIV1 that only monotonic global 
stability can be established. Using the EZ method, we show 
in section [V] that there is a region of parameter space within 
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which Kolmogorov flow is globally stable, but not monoton- 
ically globally stable. The solid and dashed curves in Fig. [T] 
illustrate this situation. 

A flow is said to be linearly unstable if the linearized ver- 
sion of Eq. ^ has exponentially growing eigensolutions. This 
is the subject of section Unl 



m. LINEAR INSTABILITY 

Our goal in this section is to find the "neutral surface" in 
the (/j,, 0, ^)-space below which ip\_{x) is linearly unstable. 
We linearize Eq. (O and write ip in Floquet form [12, 14]: 



(12) 



where —1/2 < fc < 1/2 and (p(x) has the same periodicity as 
ip\_(x), i.e., <f(x) — Lp(x + 2ir). The resulting eigenproblem is 



ii/(D z - I 2 ) 2 + I smx{D' z - X* + 1) - i0D 



= X(D 2 - l 2 )ip 
where the differential operator D and the eigenvalue A are 



(13) 



D = — + ifc . 

ax 



\ = 



The two parameters in Eq. (fT3l are defined as 



(3 



(14) 



(15) 



The basic flow is linearly unstable if there exists an u> with 
positive imaginary part, i.e. = a3{A} — fj, > 0. 



A. Tracing the neutral surface 

For a given pair of (0, v), we solve Eq. ([T3l l numerically to 
obtain the eigenvalue spectrum {X n (k, I; 0, vj} as a function 
of the wavenumber (fc, I). The integer n on A„ indexes the 
eigenbranch. Once we possess the spectrum, we define the 
function 



Anout(/3, v) = max3{A„(fc, I; 0, u)} . (16) 



Figure |2] shows max„ 3{A ra (fc, I; 0, i>)} at four values of \x 
and which happen to fall on the inviscid neutral curve. 
The function jl nen t(0 , v) in Eq. (IT6b is obtained by searching 
through the (fc, Z)-plane in Figure|2]to find the maximum. 

Because a// nou t — is the growth rate of the most unstable 
mode, the neutral surface in the (/x, 0, z/)-space is defined by 



Given fi. neu t(0, v), the neutral surface is traced using the para- 
metric equations: 



(18a) 
(18b) 
(18c) 
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The expressions above follow from Eq. ©, Eq. ( fT51 ) and 
Eq. ( fT7T >. The neutral surface is obtained from Eq. d 1 8ab 
through Eq. d 1 8cb by varying the parameters (0, v) between 
and oo. 



B. Gill's inequality 

In the inviscid case, v — 0, an argument of Gill and 
Lorenz II 1311 shows that it is sufficient to search for unstable 
modes within the circle, 



fc 2 +l 2 < 1 



(19) 



(J, = ajl neu t{0,v) . 



(17) 



In Appendix [A] we show that Eq. ( fT9b also applies to lin- 
early unstable modes of the viscous problem. Gill's inequal- 
ity underscores Fj0rtoft's result [ 1] by showing that a growing 
eigenmode must transfer energy to spatial scales which are 
larger than that of the laminar flow in Eq. ©. 



C. The inviscid case, v = 

The procedure outlined above is particularly simple in the 
inviscid case v = 0. The inviscid neutral curve in the (/i, 0)- 
plane is the solid curve in Fig. [3] The neutral curve intersects 
the /i-axis at \i = 0.52. After disentangling notational dif- 
ferences, /j — 0.52 agrees with Thess's value for unbounded 
Kolmogorov flow 111 211 . 

An interesting feature of the neutral curve is the appearance 
of a kink at the point K:{p, 0) w (0.37, 0.46). This is due to 
a change in character of the most unstable mode, evident be- 
tween the first and third panels of Fig. [2] As one moves along 
the neutral curve by decreasing fi, the peak eigenvalue on the 
k = line is getting smaller and is eventually overtaken by 
another peak emerging at a location with fc ^ 0. Consequently 
the wavenumber k of the most unstable mode discontinuously 
jumps from fc = to fc = 0.15 at K. Figure|4]shows the tran- 
sition by plotting the variation of /2 ncu t and the wavenumber 
of the most unstable mode along the neutral curve. 

We find that for all growing modes, except for those with 
fc = 0, the real part of uj is non-zero, 5R{cj} ^ 0. Hence 
the jump in fc to a non-zero value as n decreases also signals 
that the most unstable mode becomes a traveling wave. In 
summary, along the neutral curve in Fig. [3] the most unstable 
modes to the right of point K are stationary disturbances with 
k = 0, and those on the left of point K are traveling waves 
with fc 7^ 0. 
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13=0.19, ^1=0.49 



[3=0.47, ^=0.37 



[3=0.59, ^=0.33 



13=1.22, ^=0.19 




0.50 



FIG. 2: (Color online) This shows max„ Ss{A„(fc, I; /3,D — 0)} as a function of k and /; dark areas are regions of large values. The values of 
(fj,, (3) in the titles are obtained by substituting and the peak value, /incut, into Eq. dl8at and Eq. l !18bt . 
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FIG. 3: Results from linear stability analysis of Eq. (TJl. For the 
inviscid case v = 0, the solid line is the neutral curve. The point K 
is where the most unstable mode changes character. For v = 0.1, the 
dashed line is the neutral curve. 
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FIG. 4: (Color online) Variation of (a) /incut, defined in Eq. J 1 6b . and 
(b) wavenumber (k, I) of the most unstable mode along the linear 
stability neutral curve. The dashed lines in (a) show the piecewise 
linear approximation Eq. d21t . 



We mention that in the inviscid case, the neutral curve can 
be written as 



13= J*%L- M 2. 



(20) 



and Fig.Ua) indicates that 

Mneut(M) ~ 0.26- 0.12/i, 0</i<0.37, (21a) 
/w(m) ~ 0.09 + 0.34/i, 0.37 < ^ < 0.52. (21b) 



D. The effects of viscosity 

When v ^ 0, the neutral curve can be traced out by a sim- 
ilar procedure. The dashed curve in Fig. [3] is the result for 
v = 0.1, showing the stabilizing effect of v. For v = 10~ 3 , 
the neutral curve is indisguishable from the inviscid curve, in- 
dicating the limit ^ — ^ is non-singular. We find a kink K 
in the neutral curve, and a corresponding change in the struc- 
ture of the most unstable mode, always appears at (3 ~ 0.46, 
independent of v. 



IV. THE ENERGY METHOD 



Substituting Eq. ( |2TI ) into Eq. ( |20b gives an approximate ex- 
pression for the v — neutral curve; this approximation 
matches the numerical result to within the line width in Fig. [3] 



We now turn to the stability of the basic flow Eq. ([3]) us- 
ing the energy method. Fukuta and Murakami 11911 have pre- 
viously considered this problem with (3 = 0. It is easy to 
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adapt their results to nonzero [3 because the disturbance en- 
ergy equation does not involve (3: (3 enters the energy method 
only through the laminar amplitude, a((3, fi, v) in Eq. (@). 

We begin by rewriting the disturbance energy equation 
Eq. I® as 



dff y 
dt 



= 2 {aR[<p} -n)E v - 2uZ v , 



where 



(22) 



(23) 



The homogeneous functional R represents the transfer of en- 
ergy by Reynolds' stresses between the basic flow and the per- 
turbation flow. 

The basic idea of the energy method is that if the dissipative 
parameters fj, and v are large enough, then the right hand side 
of Eq. (1221 is negative for all possible i/j's (not just tp's which 
happen to satisfy the equations of motion). This implies that 
the disturbance energy decreases monotonically to zero, and 
that the laminar solution is monotonically globally stable. 

For instance, since 



(I^IKI) i 

m ~ (|V^P> "2 



it follows from Eq. (1221 that 



dt 



< 2 



E„ 



and Gronwall's inequality implies that 



(24) 



(25) 



(26) 



That is, if a < 2/i then the laminar flow in Eq. (f3]) is monoton- 
ically globally stable. This conclusion relies on the seemingly 
crude inequality in Eq. ( l24l and one expects that a stronger 
condition for monotonic global stability might be obtained by 
working harder and solving the variational problem suggested 
by maximizing the right hand side of Eq. d22l . In this spirit 
we recapitulate some of the variational results in IU9I1 by dis- 
cussing the simplest case v = in some detail. 



A. The inviscid case, v = 

With v = 0, let "Re be the maximum of R[f>] over all func- 
tion tp(x) satisfying the periodic boundary conditions: 



Re = maxR[ip] . 
•p 



(27) 



According to the energy-stability method, the function <pe{x) 
which maximizes R is the "most dangerous disturbance", 
meaning the most efficient energy-releasing disturbance. 

Applying Gronwall's inequality to Eq. ( 1221 . if clRe < M 
then Em decays monotonically to zero and the laminar so- 
lution is monotonically globally stable. Hence with v = 0, 



Ao 




FIG. 5: (Color online) Eigenvalue of the gravest mode of Eq. d29t , 
Ao(fe, /). The minimum is at k — and I — » oo. The inset shows a 
cut along k — (solid line); Ao(0, 1) decreases monotonically from 
2-\/2 at I — to 2 as I — » 00. The dashed line is the approximation 
Eq. _ 



the energy-stability boundary in the (fj,, /3)-plane is given by 

jj, = clRe, or 




(28) 



Since R[<p] is a homogeneous functional, we can find Re by 
maximizing ((fixtfiy cosx) subject to the constraint ^|V</?| 2 ) 
1 This leads to the Euler-Largrange equation, 

(D 2 - l 2 )<f = i/A (^cosxD - isina^ip, (29) 

where D is the differential operator defined in Eq. ( fT4l and the 
Lagrange multiplier appears as a real eigenvalue A. Then Re 
is given by 1 / A m ; n , where A m i n is the minimum eigenvalue. 

Again we represent (p using the Floquet form in Eq. ( fT2l 
(but now with u> = 0) and solve Eq. ( 1291 numerically to obtain 
Ao(fc, I), the eigenvalue of the gravest mode, as a function of 
k and /. Figure |5] displays Xo(k, I) above the (k, ^)-plane and 
the conclusion is that 



A„ 



min Anffc, I) 

k.l 



(30) 



with the minimum achieved at k — and I — > 00. The numer- 
ical result Eq. ( T30b is supported by the analysis of Eq. d29l in 
AppendixlBl 

The conclusion Re = 1/A m j n = 1/2 is anticipated pre- 
cisely by the simple inequality in Eq. d24l : in the inviscid 
case the variational solution does not improve the energy sta- 
bility boundary at all. The only reward from the variational 
solution is that it provides the form of the most dangerous dis- 
turbance tp E {x). The analysis in Appendix 151 suggests that 
this disturbance can be approximated by the trial function 



<Pe(x) ~ lim cos [l(y + sin a;)] cxp ( — cos2x 
/ — >oo V 2 j 



(31) 
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FIG. 6: The functional 3? (Eq. J23t ) evaluated at the trial function 
Eq. J3U with different /. The inset shows the trial function with 
I = 4. 

Fig.[6]shows that 31 evaluated at the trial function above indeed 
approaches one-half as I — > oo. 

The Reynolds' stresses of the most dangerous disturbance 
are concentrated in the neighbourhood of cos x = ±1. Since 
there is no penalty attached to using very small spatial scales 
in the functional 9l[(p], this concentration can be made ever 
more intense by taking I — > oo. The trial function shown 
in the inset of Fig. [6] illustrates this strategy for maximizing 
Jl[(p]. But <pe(x) is realized only in the I — > oo limit, and 
this is distinctly unphysical. This point is discussed further in 
Section[V] 

With Ji E = 1/2, it follows from Eq. (|28]) that the inviscid 
energy-stability boundary is given by 

0= \j^-» 2 - (32) 

Eq. ( 1321 ) is plotted as a dotted curve in Fig. [7] Despite the 
pathology discussed above, it is remarkable that the energy- 
stability boundary parallels the linear-stability neutral curve 
of Section [III] even though the parameter /3 enters the energy 
method only through the function a((3, /Lt, v) in Eq. (gj). This 
is an indication that the grossest features of the instability are 
determined by the amplitude of the laminar solution and the 
disturbance energy equation Eq. (1221 . 

B. The viscous case 

Turning now to v ^ 0, the main point of interest is how 
small viscosity affects the energy-stability curve in Fig. [7] In 
this case the solution of the variational problem can provide a 
substantial improvement over the simple inequality Eq. d24b . 
Following Fukuta and Murakami 1 19], the viscous version of 
Eq. d29l can be solved numerically and the results are shown 
in Fig. [8] As v — > the energy-stability curves limit to the in- 
viscid curve in Eq. d32l . Thus the limit ^ — ^ is not singular. 




FIG. 7: A summary showing three different stability boundaries for 
two-dimensional Kolmogorov flow. The dashed curve is the neutral 
curve of the linear stability documented in Section 11111 the dotted 
curve is the energy stability boundary in Eq. d32fc ; the solid curve 
is the EZ stability boundary in Eq. d40| >. The point K at (n, (3) = 
(0.37, 0.46) is where the most unstable linear mode changes from a 
stationary disturbance to a traveling wave. 




FIG. 8: This shows energy-stability boundaries in the (fi, (S) parame- 
ter plane. The flow is monotonically globally stable if the parameters 
(/3, n, v) locate it above the appropriate curve. 

One other point of interest is the wavenumber I* of the most 
dangerous disturbance. Fig.|9]shows the dependence of Z* on 
v. In agreement with the small- wavenumber expansion lfl9Tl . 
if u/fi > 0.0204 then U = and He = 1/^8(1 + ^). 
Additional asymptotic analysis of ours (not presented here) 
shows that if u/jj, < 0.0204 then 

l * = (i) 1/3+0 ^- (33) 
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FIG. 9: This shows the dependence of the wavenumber Z* of the 
most dangerous disturbance on v. (The most dangerous disturbance 
has k, = for all v.) If v > 0.0204/1 then U = 0. As v decreases, 
U increases and ultimately, as show in the insert, h ~ v~ 



In the limit v — > we recover the inviscid result U 



V. THE ENERGY-ENSTROPHY METHOD 

The gap between the neutral curve and the energy-stability 
boundary in Fig. [7] reflects the stark difference between the 
linearly unstable eigenmodes and the most efficient energy- 
releasing disturbance identified by the energy method. Specif- 
ically, the wavenumbers of the exponentially growing linear 
modes satisfy Gill's inequality Eq. ( fT9] l while, according to 
the energy method, the most efficient energy-releasing distur- 
bance has a wavenumber (k, I) = (0, oo). To improve on the 
energy method, the energy-enstrophy (EZ) method uses both 
the disturbance energy and enstrophy. 



A. The EZ method, v = 

With v = 0, Eq. (TO) implies that 

E v (t) - Z v (t) = e -2 "* [E V (Q) - Z v (0)} . (34) 

Thus, although an arbitrary initial condition may not fall in 
the set 



EZ 



x) : E ir _ 



(35) 



the solution to Eq. (|6]l (with v=0) is attracted by EZ. And if 
the initial disturbance does happen to fall within EZ then it 
stays within EZ. We refer to initial conditions which fall in 
EZ as ".EZ-disturbances". 

Figure [10] illustrates the approach to the attracting set EZ 
using a numerical solution of Eq. ©. Panel (a) shows that 



solutions starting outside of EZ are quickly attracted to EZ; 
panels (b) and (c) show that solutions starting in EZ, remain 
close to EZ. The trajectories with fi = 0.4 in panels (a) 
and (b) show two unstable solutions; nonlinearity halts the 
growth of the disturbance so that both trajectories asymptote 
to (E V .Z V ) w (0.49,0.49) as t — > oo. The solutions with 
fi = 0.61 in panels (a) and (c) show disturbances whose en- 
ergy decays monotonically to zero (although \i = 0.61 is be- 
low the energy stability boundary). These numerical results, 
with < v -C /J, show that EZ is an attracting set with small 
but nonzero viscosity. 

To obtain an energy-stability bound for ^^-disturbances 
we seek 



JIez = max !R\cp] . 



(36) 



This is equivalent to maximizing 0i[<p] in Eq. d23l subject to 
the constraints 



flvyf) = <(W) a > = i 



(37) 



The solution ifiEz(x) is the most dangerous UZ-disturbance. 

The search in Eq. d36] l is over a smaller set than the one used 
to define He in Eq. d27b , and so 31ez < $-e- In anticipation of 
the variational calculation in the next subsection, we remark 
that 

%ez - 0.3571 <R E = \. (38) 

Furthermore, the argument surrounding Eq. ( IA7I ) shows that 
the Floquet wavenumber of an £JZ-disturbance must satisfy 
Gill's inequality Eq. f!9) . Thus the pathology of the energy 
method is avoided. 

The disturbance energy of an SZ-disturbance satisfies 



dEf 
At 



(39) 



and Gronwall assures us that the energy decays monotonically 
to zero provided that a3? ez < M- Using the definition of a in 
Eq. ©, this condition leads to the EZ stability curve 



'3* _ ^2 



(40) 



shown as a solid curve in Fig. [7] The trajectories with // = 
0.61 in panels (a) and (c) of Fig. [TOl show two solutions just to 
the right of the EZ stability curve in Fig. [7] 

In Fig. [7] there is a gap between the EZ curve and the 
energy-stability curve. In this region SZ-disturbances de- 
cay monotonically, while general initial disturbances might 
have amplifying disturbance energy. In fact, initial amplifi- 
cation can be assured by deploying the trial function on the 
right of Eq. OTb at t = 0. We argue heuristically that this 
amplification can only be transient: Eq. d34l l shows that the 
set in Eq. fl35l l attracts all initial conditions. Thus, since EZ- 
disturbances decay monotonically to zero, general initial dis- 
turbances will also eventually decay to zero as they evolve 
towards EZ via Eq. d34l l. We are arguing heuristically that in 
the gap between the EZ curve and the energy-stability curve 
the laminar solution is globally stable, but not monotonically 
globally stable. 
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FIG. 10: (Color online) Trajectories of some illustrative solutions of 
Eq. (6) in the \E v (t), Z^(t)]-plane; all solutions have /3 = and 
10~ 3 < v < 5 x 10~ 3 . The dashed line marks the set EZ in 
Eq. (35). 



B. The EZ variational problem 

To obtain Hez = 0.3571 in Eq. ([38), we maximize 
(ip x ip y cos x) subject to the constraints that the disturbance 
enstrophy and energy are both equal to one-half. A direct 
assault based on solving the Euler-Lagrange equation is de- 
scribed in Appendix ICl It is more instructive to develop some 
intuition using the trial function 



for which 



ip = Aq cos ly — Bi sin ly sin x . 



((p x ip y cosx) = -lA Bi . 



(41) 



(42) 



With the constraints in Eq. ( 137) . the parameters Aq and B\ 
are expressed in terms of the wavenumber / and substituted 
into Eq. (02). Thus the trial-function in Eq. fiD leads to the 
estimate 



Hez ~ max I 



(1-/2) 



I \2{\ + l 2 ) 



(43) 



The maximum of the right hand side is 1 — 2 1 I 2 = 0.293 • • • , 
which is achieved at I 2 — \/2 — 1. This is a lower bound on 

To more closely approximate Hez we generalize Eq. (HP 

to 

oo oo 

ip = cos ly A 2n cos 2nx— sin ly B 2n +i sin(2n + l)a; . 

n=0 n=0 

(44) 

TABLE J] summarizes the result of retaining more terms in 
the Fourier series Eq. d44) . Optimization over A n , B n and I 



FIG. 11: The most dangerous i5Z-disturbance Eq. (44) calculated 
using the final row of TABLE|I] 



A 2 



A 4 



Bi 



B, 



I 



■31 



EZ 



1.4142 — 
2.1561 0.2764 
2.2141 0.2694 



1.2872 — 0.6436 0.2929 
1.3010 — 0.4253 0.3556 
1.2983 0.0371 0.4171 0.3571 



2.2175 0.2689 0.0037 1.2984 0.0376 0.4166 0.3571 

TABLE I: The estimate of 3?bz in the final column quickly con- 
verges to the numerical value indicated in Eq. d38b as more terms in 
the Fourier series Eq. (44) are retained. 



provides an estimate of OIez that quickly converges to the nu- 
merical value indicated in Eq. (1381 . Fig.[TT]illustrates the most 
dangerous SZ-disturbance calculated using the final row of 
TABLE U 



VI. CONCLUSION AND DISCUSSION 

We have studied the instability of geophysical (f3 ^ 0) Kol- 
mogorov flow and focussed on the limit in which drag is much 
stronger than viscosity, fi ^ v. We show that the form of the 
fastest exponentially growing eigenmodes changes abruptly 
from stationary disturbances to traveling waves as one moves 
along the neutral curve in the (/x, j3) plane. In this sense, 
P strongly affects the instability. Nonetheless, the stability 
boundary obtained by the energy method lfl9tl is reasonably 
close, and parallel, to the linear-stability neutral curve. Be- 
cause (3 enters the energy method only through the amplitude 
of the laminar solution i.e., et(/3, /i, v) in Eq. (0), this success 
indicates that the main features of the instability are deter- 
mined by a and the disturbance energy equation. 

It is striking that the most unstable inviscid disturbance 
identified by the energy method is unphysical because it is 
strictly realized only in the infinite wavenumber limit, I — > oo. 
Despite this issue, the energy method still delivers a useful 
sufficient condition for global monotonic stability. 



9 



In section [V] we extend the energy method by incorporat- 
ing information based on the disturbance enstrophy equation. 
This EZ method results in a third stability boundary, the EZ 
boundary, that lies closer to the neutral curve than does the 
energy-stability boundary. One might therefore conclude that 
the additional information provided by enstrophy produces 
only a quantitative narrowing of the gap between energy sta- 
bility and linear stability. However the main interest in the EZ 
method is that it identifies a particular type of slowly evolving 
disturbances i.e., those disturbances belonging to the set EZ 
in Eq. (|35t . The EZ method relies on identifying the most 
unstable disturbance in EZ, and this disturbance, given ap- 
proximately by the Fourier coefficients in TABLEU has finite 
enstrophy. 

The various types of stability we have discussed here char- 
acterize the neighbourhood of the laminar solution ip\_ (x) in 
Eq. OJ. However Doering and Constantin ll25l l26ll have re- 
cently devised a variational procedure in which notions of 
energy stability are applied to statistically steady turbulent 
flow. This results in bounds on important large-scale quan- 
tities, such as the mechanical energy dissipation and the heat 
flux. A motivation for this paper has been the possibility of 
applying or generalizing the technique of Doering and Con- 
stantin to two-dimensional turbulence. In this context it is es- 
sential to take account of enstrophy conservation, perhaps via 
the notion of EZ-stability. 
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The first factor is non-negative, and consequently a growing 
eigenmode must have more energy than enstrophy: 



E v (0) > Z v (0). 



(A5) 



In the inviscid case considered by Gill lfl4ll . the result Eq. ( IA5b 
is an equality following immediately from Eq. (lAlb without 
the Cauchy-Schwarz excursion. 

To obtain the inequality Eq. ( [T9| >, substitute the Floquet 
form Eq. (fT2l . with 



into Eq. iA5i : 



y{x)= J2 Vne m \ (A6) 



J2 (1-/4KKI 2 >0. (A7) 



Above K 2 n = (k + n) 2 + l 2 is the wavenumberof the n'th wave 
in the Floquet series. Recalling that —1/2 < k < 1/2, it is 
easy to see that 1 — /i 2 is negative if \n\ > 2. Consequently 
the inequality above can only be satisfied because 1 — h 2 , or 
1 — h\ x are positive. This requirement implies Eq. ( fT9] l. 



APPENDIX B: THE INVISCID ENERGY STABILITY 
EIGENPROBLEM 

In this appendix we use perturbation theory to analyze the 
eigenproblem Eq. d29l ) in various limits. The analysis can be 
simplified by the transformation fl9ll 



APPENDIX A: EXTENSION OF GILL'S INEQUALITY TO 
INCLUDE VISCOSITY 

Gill's inequality in Eq. ( [T9T > restricts the Floquet wavenum- 
ber of an exponentially growing eigenmode to lie within the 
unit circle in the (fc, Z)-plane. To prove this result, observe 
that a growing eigenmode of the Floquet form Eq. (fT2l satis- 
fies Eq. ([Toll. Thus 

(w, + n) [E V {<S) - Z v (0)} = isP v (0) - vZ v {Q) , (Al) 

where u>.- L > is the imaginary part of u> and E v (0) et cetera 
is the initial energy et cetera of the eigenmode. 
The disturbance enstrophy can be written as 



so that the Cauchy-Schwarz inequality implies 

Z 2 < P^Eu, . 



(A2) 



(A3) 



Replacing P v (0) in Eq. (TAJ) by Z 2 (0)/E v (0) we obtain an 
inequality equivalent to 

{ LJl + ^ l + U E~§j) [^(°)--M0)]>0. (A4) 



tp(x) — exp 



i^Asin x 



0(x) 



This results in 

— 2ifc— 
dx 2 dx 



cos x — 1 



k'< 



(Bl) 



= 0. (B2) 



Fukuta and Murakami lfl9ll solve Eq. (1B2I) for the gravest 
eigenmode in the small wavenumber limit: 

9 = 1 + -(fc 2 + 1 2 ) cos 2x - iifc(fc 2 + I 2 ) sin 2x + 0(k, I) 4 , 

4 " (B3) 

with the corresponding gravest eigenvalue 

This shows that as (k, I) -> (0, 0), with I ^ 0, A -> 2^2. 

We now consider the complementary case k — and I 3> 1; 
we define a small parameter by e 2 = Numerical solu- 
tion of Eq. (IB2b indicates that the minimum of Xo(k, I) is at 
Ao(0,oo) = 2 and we are seeking some analytic assurance 
of this hypothesis. The eigenfunction is concentrated in the 
neighbourhood of x — where 4- cos 2 x — 1 in Eq. iB2\ is 
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slightly positive. Thus we expand cos 2 x ~ 1 — x 2 + 0(x A ) 
and introduce a boundary-layer coordinate X — x/e so that 
Eq. (IB2I) reduces to 



'XX 



Ae 2 



-X 



>\e = o (e 2 ) . 



(B5) 



We recognize the quantum harmonic oscillator equation and 
thus obtain the eigenvalues as 



A 2 -4 



= 2n 



n = 0,1,2, 



(B6) 



2A„e 2 

Hence, the eigenvalue of the gravest mode (n — 0) is 

A «2 + r 1 + O(r 2 ) , (B7) 

and Ao — * 2 as I — ► oo. 
The expression 



\2 



\/l + c 2 l 2 
1 + cl 2 



= 1 - 2~ 1/2 = 0.2928- •• 



(B8) 

agrees with Eq. ( IB4I ) when k = and I <C 1, and with 
Eq. (IB7t if k = and I ^> 1. The inset of Figure [5] shows 
that the interpolation Eq. dB8b is a good approximation to the 
numerically computed eigenvalues for all I. 



APPENDIX C: THE EZ EULER-LAGRANGE EQUATION 

The direct approach to solve the EZ variational problem in 
section [VB] is to include the two constraints Eq. (|37T > using 
two Lagrange multipliers p and q. Thus setting the variational 
derivative of 

(<p x <p y cosx) - q <|V^| 2 > +p [<|V<^| 2 - (V 2 ^) 2 )] (CI) 



to zero results in the Euler-Lagrange equation 

pV A <p + (p — q)y 2 f = — (f y s'mx — (p xy cosx . (C2) 

Multiplying Eq. (IC2t by ip and taking the spatial average, we 
deduce that 



(C3) 



Writing ip in the Floquet form Eq. ( TT2l (with oj = 0), we solve 
Eq. dC2b by regarding p as an eigenvalue and q as a parameter, 
for different values of k and I satisfying the Gill's inequality 
Eq. ([T9T >. For a given q, admissible solutions are the subset of 
eigenfunctions of Eq. iC2\ that satisfies 



(IV^I 2 ) 
<(W) 2 ) 



1 = 0. 



(C4) 



JIez, defined in Eq. d38i l. is then given by the maximum q at 
which such solutions exist. Now, we know that 

0.3571 < q < 0.5. (C5) 
The upper limit is what we get if we ditch the enstrophy con- 
straint and perform the maximization within the larger class of 
functions satisfying only (|V(^| 2 ) = 1. The lower limit is ob- 
tained from the trial function method described in section lVBl 
Thus, we only need to search for 'Rez within the range in 
Eq. 0). The result is 5i EZ = 0.3571 with p = 0.0069, this 
occurs at k = and I = 0.4166. These numbers are virtually 
the same as the trial function result given in the final row of 
TABLE U We observe that in both the energy method and the 
EZ method, the most unstable disturbance has wavenumber 
k = 0, we have not been able to prove this result analytically. 
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